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Abstract 

We study the chiral phase transition of QCD at finite temperature and den- 
sity by numerically solving Schwinger-Dyson equation for the quark propa- 
gator with the improved ladder approximation in the Landau gauge. Using 
the solution we calculate a pion decay constant from a generalized version of 
Pagels-Stokar formula. Chiral phase transition point is determined by analyz- 
ing an effective potential for the quark propagator. We find solutions for which 
chiral symmetry is broken while the value of the effective potential is larger 
than that for the chiral symmetric vacuum. These solutions correspond to 
meta-stable states, and the chiral symmetric vacuum is energetically favored. 
We present a phase diagram on the general temperature-chemical potential 
plane, and show that phase transitions are of first order in wide range. 
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I. INTRODUCTION 



In QCD an approximate chiral symmetry exists at the Lagrangian level and it is sponta- 
neously broken by the strong gauge interaction. Accordingly approximate Nambu-Goldstone 
(NG) bosons appear, and pions are regarded as the NG bosons. In hot and/or dense mat- 
ter, however, quark condensates melt at some critical point, and the chiral symmetry is 
restored (for recent reviews, see, e.g., Ref. |JJ). To study the chiral phase transition we need 
non-perturbative treatment. 

Schwinger-Dyson Equation (SDE) is a powerful tool for studying the qualitative structure 
of chiral symmetry breaking of QCD. By suitable choices of the running coupling, in which 
the asymptotic freedom is incorporated, the SDE provides a non-trivial solution for the mass 
function. The behavior of the mass function at high energy is consistent with that given 
by the operator product expansion technique (for a review, see, e.g., Ref. 0). The SDE is 
understood as a stationary condition of an effective potential for the quark propagator, and 
the chiral broken solution is energetically favored by the effective potential 

The SDE was applied for studying the chiral phase transition of QCD at finite tempera- 
ture and/or density |f7HIB[|. In many of those the chiral phase transition point was determined 
from the point where the mass function (or pion decay constant) vanished without studying 



an effective potential. In Ref. [jn| an effective potential for the mass function was analyzed 
by using an approximate form of momentum dependence of the mass function. It was shown 
that there were meta-stable states where the chiral symmetry was broken while the energy 
was bigger than that for the symmetric vacuum. It is interesting to study whether this kind 
of meta-stable state is derived by fully solving the SDE. 

In this paper we study the chiral phase transition in QCD at finite temperature and/or 
density by numerically solving the SDE using Matsubara formalism. As an order parameter 



we calculate a pion decay constant using a generalized version of Pagels-Stokar formula [19 



We calculate the value of the effective potential at the solution of the SDE to determine the 
chiral phase transition points. We show that meta-stable states actually exist in wide range. 
The order of phase transition are determined by comparing the order parameter with the 
effective potential. The phase transition is of the first order for wide range at T ^ and 

This paper is organized as follows. An effective potential for the quark propagator is 
introduced at finite temperature and/or density in section |1[ The SDE is derived as a 
stationary condition of the potential. In section |T| a formula for calculating a pion decay 
constant is presented. Section [IV] is a main part of this paper, where numerical solutions 
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of the SDE are shown. Values of the pion decay constant and the effective potential are 
calculated from the solutions, and the phase structure is studied. Finally a summary and 
discussion are given in section M. 



II. EFFECTIVE POTENTIAL AND SCHWINGER-DYSON EQUATION 

At zero temperature the effective action for the quark propagator S F is given by 

T[S F ] = iTr Ln S F - Tr [i?S F ] - m 2PI [S F ] , (2.1) 

where k 2 pi stands for the two particle irreducible (with respect to quark-line) diagram con- 
tribution. In this paper we take the first order approximation for K2PI, in which one gluon 
exchange graph contribute, 

/V /V C c 

«2pi[Sf] = i ~^ r 1 J d 4 xd A yghi [S F (x - y)i^S F (y - x)i lv \ D^(x - y) , (2.2) 

where the number of flavor (iVf) and of color (Ac) and the second Casimir of SU(A*c) (C2) 
appear since we factor color and flavor indices from the quark and gluon propagators. In 
this expression, the quark propagator Sp carries bispinor indices and the gluon propagator 
D^ v carries Lorenz indices. From the above effective action the effective potential for S F is 
obtained by the usual procedure. In momentum space it is expressed as 



j4 

V[3 f ] = J -r^{\ndet [S F (p)]+tr \jfiS F (p)]} 



l(27T) 

+ 1 J J^y 4 J -§^°29 2 (P, k)tr [zS F (ph^Sp(k) lu ] iD^(k - p) , (2.3) 

where an overall factor NqN{ is dropped. The running coupling g 2 (p, k) is introduced to im- 
prove Schwinger-Dyson Equation (SDE) to include aymptotic freedom of QCD. An explicit 
form of the running coupling will be given below. In this paper we take the Landau gauge 
for the gluon propagator: 



iD^(l) = i 



r- — 



(2.4) 



P 

General form of the full quark propagator is expressed as 

iS F \p) = A{p)p / -B{p) . (2.5) 

SDE for the quark propagator Sf is given as a stationary condition of the above effective 
potential j3J (See also, Ref. ||). The SDE gives coupled equations for A and B. If the 
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running coupling g(p, k) is a function in p 2 and k 2 , i.e., it does not depend on p ■ k, it is 
easy to show that A(p 2 ) = 1 is a solution of the SDE in the Landau gauge with the ladder 
approximation. 



Let us go to non-zero temperature and density. In the imaginary time formalism [£0] the 



partition function is calculated by the action given by (See, for reviews, e.g., Refs. pl| , |22| .) 

-l/T 



s 



dr Id x 



QCD 



+ /i^7 i> 



(2.6) 



where T and /i denote temperature and chemical potential associated with quark (baryon) 
number density, and £qcd takes the same form as QCD Lagrangian at T = /i = 0. From 
the above action we obtain an effective potential for the quark propagator S F similar to the 
one in Eq. 

d 3 p 



V[S f 



u 

■ d 

(2 



{lndet [S>(p)]+tr tyiS F (p)}} 



(2tt) 3 

+\ t2 £ / ^ / (0^ 2 (P, fc) tr [iS F (p)j»iS F (k)j„] uy{k - p) 



(2.7) 



where 



Po = iu + p, , -u = (2m + l)7rT 



k = iv + fi 



(2n + l)7rT , (n,m = integer) 



(2.8) 



and J2 U an< i imply summations over m and n, respectively. 

Argument of the running coupling should be taken as (k — p) 2 for preserving the chiral 
symmetry 



However, as is shown in Ref. ||23|| , angular average, i.e., the running coupling 
is a function in — p\ + |p| 2 — k\ + |fc| 2 , is a good approximation at T = [i = 0. If we naively 
extend this approximation, the running coupling will depend on \i. However, the running 
coupling should not depend on /i. Then in this paper we use the approximation where the 
running coupling is a function in —(po — /i) 2 + \p\ 2 — {ko — fi) 2 + \k\ 2 . The explicit form of 
the running coupling is pi 



g 2 (p,k) = (47rY 



lliVc - 2N f 



A In 



t 



' u 2 + x 2 + v 2 + ?/ 2, 
if i F < t 



^ tp ^tp 



(2.9) 
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with 



x = \p\, y=\k\. (2.10) 

A qc d is a scale where the one-loop running coupling diverges, and the value of A qc d will be 
determined later from the infrared structure of the present analysis.^ tp and tc are parame- 
ters introduced to regularize the infrared behavior of the running coupling. In the numerical 
analysis below, since the dominant part of the mass function lies below the threshold of 
charm quark, we will take Nf — 3 and iVc = 3. Moreover, we use fixed values tp = 0.5 and 
t c = 2.0 [§. 

General form of the full quark propagator, which is invariant under the Parity- 
transformation, is given by [f7],|ll||12l,|15| 



iS^(p) = A(p)j- B(p) + C(p) Pol ° - % - [ 7 ° , p] D{p) . (2.11) 

Here A ~ D are functions in p and |p| (as well as T and /i), although we write A(p), etc. 

It should be noticed that the effective potential in Eq. (|2.7| ) is invariant under the tran- 
fromation: 

5f(po,p) — S' F {po,p) = (il 1 ! 3 ! ) SUPo, -P) (^V) , (2-12) 

where Sp implies transposition of in spinor space.0 At T — \i — this is nothing but 
the time reversal symmetry. It is natural to expect that this "time reversal" symmetry is 
not spontaneously broken. Existence of this "time reversal" leads D(p) =0. 

Moreover, at fi = the effective potential is invariant under the following "charge con- 
jugation" : 

S F (p) — S' F (P) = -(*7 2 7 )^(-p)(*7 2 7 ) • (2.13) 

We also assume this symmetry is not spontaneously broken, and it implies that all the scalar 
functions A, B and C (as well as D) are even functions in p = iu. As we mentioned above, 
at T — [A — 0, of course C{p) = 0, and A{p) = 1 in the Landau gauge. C{p) ^ as 
well as A(p) ^ 1 does not imply chiral symmetry breaking, and non-zero B(p) is the only 
signal of the breaking in the present analysis. Then we consider A{p) — 1 = C(p) = as 



1 The usual Aqcd is determined from the ultraviolet structure. 

2 In our convention, (7 )* = 7 , (7 1 )* = 7 1 and (7 3 )* = 7 3 , while (7 2 )* = — 7 2 . 
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an approximate solution at general T and //. We note that non-zero D(p) also breaks chiral 
symmetry, however this term do not contribute to local order parameter of ciral symmetry, 
(ipip) = — J trSp- We choose vacuum of "time reversal" invariance as discussed above.f] 

As in the case at T = /i = 0, the SDE is given as a stationary condition of the effective 
potential (|2"7I 



f d k 

i^(p)-y=rj] / -— 3 C 2 g 2 (p,k) lfli S F (k f ) lui D^(k-p) . (2.14) 

V (27T) 

By taking a trace and performing the three-dimensional angle integration (Note that the 
present form of the running coupling does not depend on the angle.), we obtain a self- 
consistent equation for B: 

where 

^.^.^( fezglglg ). ,, 16) 

In general cases B is a complex function. By using the fact that K(p, k) is a real 
function, it is easily shown that B*(p) satisfies the same equation as B(p*) does, where 
p*fj. _ (jP* ^ = (_2-u + So these two are equal to each other up to sign, B*(p) = B(p*) 
or B*(p) = —B(p*). Since B is an even function in p at \x = 0, the choice of positive 
(negative) sign implies that B is real (pure imaginary) at /i = 0. B should be real at /x = 0, 
then 

5» = 5(P*) • (2.17) 

After a solution of the SDE fl2.14| ) is substituted with the approximate solution A(p) — 1 = 
C(p) = into Eq. (|2.7p , the effective potential becomes 



V[B sol ]=V[B sol ]-V[B = 0] 



71 



In solVf " r + 



-P 2 / 5 soi(p) - P 



(2.18) 



where 5 so i denotes a solution of Eq. ( p.!5| ). This value of the effective potential is understood 



as the energy density of the solution. So the true vacuum should be determined by studying 



3 Acctually, in the numerical calcuration of SDE in bifercation apploximation, no D(p) ^ solution 
is obtained 
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the value of the effective potential. When the value of V"[-B so i] for a nontrivial solution B so \ 
is negative, the chiral broken vacuum is energetically favored. Positive V[_B so i], however, 
implies that the chiral symmetric vacuum is the true vacuum. 



III. PION DECAY CONSTANT 

At zero temperature and zero chemical potential the pion decay constant is defined by 
the matrix element of the axial-vector current between the vacuum and one-pion state: 

(01^(0)1^))= (3.1) 

where a and b denote iso-spin indices and Jg„ = ip^^^T a il) is an axial- vector current with T a 
being a generator of SU(iVf). At finite temperature and density there are two distinct pion 
decay constant according to space-time symmetry of matirix element, which are defined 



by 

(0| J^(0)\Aq))t,, = *5 ab K(V ■ q)U L + (g, u - ^K)?"At] , (3-2) 

where = (1, 0) is a vector of medium. In this expression f n i and f n T are defined in the 
zero momentum limit, q —>■ | ]25|| . 

The above pion decay constants are expressed by using the amputated pion BS amplitude 

X as 

V,(V ■ q)U + - V,V u )q v U T 

= "§^ T E / 7^ tr [wiS F (p + q/2)x(j>;q)iS F (p-q/2)\ . (3.3) 

By using the current conservation, it is shown that the pion momentum q and the pion 
energy q Q satisfy the dispersion relation |25| 



ql = f -f \qf ■ (3-4) 

It is straightforward to obtain the pion decay constant if we have explicit expresstion of 
the amputated pion BS amplitude. Although it is generally quite difficult to obtain the BS 
amplitude, one can determine it in q — > limit, x(p; q — 0), from the chiral Ward-Takahashi 
identity 

iq^ (p - q/2, p + q/2) = Sp 1 (p - q/2) T a7s + T^Sp 1 (p + q/2) , (3.5) 

where is the axial vector-quark-antiquark vertex function, and we have suppressed color 
indices. (Note that q° in Eq. (|3.5|) is independent of q, i.e., they do not generally satisfy the 
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pion on-shell dispersion relation ( |3.4|) .) In the zero momentum limit, q, q° — > (the on-shell 
limit of the pion), is dominated by the pion-exchange contribution: 

1 



iT^(p-q/2,p + q/2) 



K(V ■ q)U + (g, u - V.VvWUt] 



ql-^l 



x(p,0)T a , (3.6) 



where u 2 = (fnr/ I-kl) \q\ 2 - Substituting Eq. ( |3.6|) into Eq. ( |3.5| ) and taking the zero mo- 



mentum limit, we find 



x(p; o) = 7— .b(p)75 

JnL 



(3.7) 



where we have used D(p) = in Sp(p)- 

The approximation adopted by Pagels and Stokar [TjJ is essentially neglecting derivatives 
of x i n the zero momentum limit (See, for reviews, e.g., Refs. [PPB"|): 



d d 
q->od\q\ g-^o dq 



(3.8) 



This approximation is same as replacing x(p] q) in Eq- (|3.3|) with x(p] 0), and reproduces the 
exact value of f n at T = in the ladder approximation within a small error 0. Then we 
use the same approximation at finite T and p in this paper. Replacing x(p; q) with x(p; 0) 
in Eq. (|3.3j ), differentiating Eq. ( [3. 3D with respect to g Q and taking g — > limit, we obtain 



—TV" 
2 ^ 



(2tt) j 



tr 



7,75 U%^*(p; o)Ufr(p) + <&(p)x(p; o)^ (p) 



2 dp c 



Here we have incorporated the on-shell dispersion relation 
into Eq. ( |3.9| ) and taking V^V a part, we find 



2 <9p a 
Then substituting Eq. 



(3.9) 



fl L = AN c TY. 



dxx 2 B{p) (B(p)-p ^ 



dp 



2tt 2 



(3.10) 



(b 2 ( p ) -p 2 y 

where we performed the three-dimensional angle integration. If the mass function B(p) is a 
function in p 2 = Pq — p 2 , this agrees with the formula derived in Ref. |TJ]]. We also obtain a 
similar formula for f nT by selecting the (g^ a — V^V a ) part: 



UtUl = in c tJ2 



dxx 2 B(p) (B(p) - 



x dB(p) 
3 dx 



2tt 2 



(B 2 (p)-p 2 Y 

where f n L in the left hand side appears from the normalization of x ( see Eq. ( |3 . 7|) ) . 
It is convenient to define a space-time averaged pion decay constant by 



(3.11) 
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V,V u )q»U T ) 



7 {flL + %UtUl 



(3.12) 



This expression aggrees with Pagels-Stokar formula jl9[ at T = and /i = after replacing 
the integral variables. 

In the formula ( |3.10| ) we need a derivative of B(p) with respect to po. However what is 
obtained by solving the SDE ( |2.15| ) is a function in discrete u = (2n + l)nT. To obtain the 
derivative an "analytic continuation" in u from discrete variable to continuous variable is 
needed. This "analytic continuation" is done by using the SDE ( |2.15|) itself: 

dB(p) dK(p,k) B(k) 



Po- 



dp 



dK(p, k) 

Po — ~ * 

opo 



B 2 {k) - k 2 ' 



(3.13) 



IV. NUMERICAL ANALYSIS 

In this section we will numerically solve the SDE (|2.15|) and calculate the effective po- 



tential and the pion decay constant f n (T) defined in Eq. (|3.12|) for the cases 1) T = and 
/i 7^ 0, 2) T 7^ and /i = and 3) T ^ and ji ^ 0, separately. The essential parameters 
in the present analysis are A qc( j and tp in the running coupling (|2.9| ). In this paper we fix 
tp = 0.5 and the value of A qc d is determined by calculating the pion decay constant through 
the Pagels-Stokar formula |B1 at T — jj, — 0. Most results are presented by the ratio to 



the pion decay constant at T — n — 0. When we present actual values, we use the value 



of the pion decay constant in the chiral limit, /„- = 88MeV [27]. (This value of /„- yields 
A qcd ~ 582 Me V.) 

A. Preliminary 

Here we will summarize the framework of the numerical analysis done below. Variables 
x and y in Eq. ( p. 10 ) take continuous values from to oo. To solve the SDE numerically 



we first introduce new variables X and Y as X = ln(x/A qc d) and Y = ln(?//A qc( i). These 
variables take values from — oo to oo. Dominant contributions to the decay constant and 
the effective potential lie around 0, i.e., x or y is around A qc( j, as shown later. Then we 
introduce ultraviolet and infrared cutoffs, X, Y e [Air, Auv]- Finally, we descretize X and 
Y at Nx points evenly: 
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Xi,Yi = A m + AX-t , i = 0,l,...,N x -l , 

AX = Auv ~ AlR . (4.1) 
N x -l V } 

Accordingly integrations over x and y are replaced with appropriate summations 

Jdx, J dy -> AX e Xl , Ax E e ^ ■ ( 4 - 2 ) 

i i 

We extract the derivative of B(p) with respect to x, which is used in the formula for f^T in 
Eq. Q3.ll ), from three points by using a numerical differentiation formula: 



dB(x) 

x 



dx 



AB{x l+l ) - 3B( Xi ) - B(x i+2 ) 



/a x- 2AX 



In the case of T = Matsubara frequency sum becomes an integration over u or v from 
— oo to oo. By using the property (|2.17| ) we can restrict this region from to oo. Then we 
will perform a procedure similar to that for x and y: 

U n = ln(w n /A qcd ) , V n = ln(w n /A qcd ) = A IR + AU ■ n , n = 0, 1, . . . , Ny - 1 , 

AU = Auv ~ AlR . (4.4) 

In the case of T ^ we truncate the infinite sum of Matsubara frequency at a finite number 

oo Njj 

E - E • (4-5) 

n=-oo n =—Nu — l 

The above Eqs. (|4.1|)-([4T5|) set regularizations adopted in this paper. We will check that the 
results are independent of the regularizations. 

Solving the SDE is done by an iteration method: 

B ncw (p) = K(p,k)* -ggggffl - . (4.6) 

Starting from a trial function, we stop the iteration if the following condition is satisfied: 



sv V ( sv 



* 



5[S F (p)]oidJ \8[S F {p)] i di 
= (B old (p) - B ncw (p))^ * (B old (p) - B ncw (p)) , (4.7) 

with suitably small e. To obtain the second line we have used the fact that A(p) — 1 = 
C(p) = is an approximate solution. This condition is natural in the sense that the 
stationary condition of the effective potential is satisfied within a required error. In this 
paper we use e = 1CT 10 . 
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B. T = and fj, ^ 



As we discussed in the previous subsection the infinite Matsubara frequency sum becomes 
an integration over continuous variable u or v at T = 0. Descretizations of these variables as 
well as x and y are done as in Eqs. ( |4.1| ) and ( |4.4| ). We use the following choice of infrared 
and ultraviolet cutoffs: 

X, Y e [-4.9, 2.9] , 

U, V e [-10.0, 2.8] . (4.8) 
For initial trial functions we use the following two types: 

(A) B(p) = B sol (p)\ T=fi=0 , 

( B ) B(p) = { 01 ' Vd ' f ° r Un ~ U2 ° ^ Xl ~ Xl ° ' (4.9) 

[ , otherwise , 

where B so i(p)\ T =Q is a solution of the SDE at T = \i = 0, and U n and Xj are descretized 
variables as in Eqs. ( |4.1| ) and ( |4.4| ). To check validity of the cutoffs in Eq. ( |4.8| ) we show (a) 
real part and (b) imaginary part of the solution in Fig. [I], and integrands of (a) /^(/x), (b) 
— V[Sgoi] in Fig. H, at fi/ f n = 3.0 for (Nu,Nx) = (70,60) by using the initial trial function 
(A). Figure |l|(a) shows that the real part becomes small above A qC( j. This behavior is similar 
to that of the solution at T = /i = 0. The dependence on x of the imaginary part is similar 
to that of the real part. Since the imaginary part is an odd function in u, it becomes zero 
in the infrared region of u. Figure [2] shows that the dominant part of each integrand lies 
within the integration range. These imply that the choice of range in Eq. (|4.8|) is enough. 

Next let us study the dependence of the results on the size of descretization. We show 
typical values of f^ip)/ in Fig. |3] and that of V[B so \]/ f% in Fig. [| for four choices of the 
size of descretization, (Nu,N x ) = (40,30), (50,40), (60,50) and (70,60). This shows that 
the choice (Nu, Nx) = (70, 60) is large enough for the present purpose. 

Now we show the resultant values of fn{/J>)/fn for two choices of initial trial functions 
in Fig. |5|. Below fi/f n = 3.0 both the trial functions converge to the same non-trivial 
solution. The value of /^(/x) increases slightly. The choice (B) in Eq. ( |4.9| ) converges to the 
trivial solution above /j,/ f n = 3.0. However, the trial function (A) converges to a non-trivial 
solution, and the resultant value of /^(/i) increases. Above \ij = 6.8 (A) also converges 
to the trivial solution. In the range 3.0 < fi/f w < 6.8 two initial trial functions converge 
to different solutions: one corresponds to the chiral broken vacuum and another to the 
symmetric vacuum. Moreover, the trivial solution is always a solution of the SDE ( |2.15 ). 
Then we have to study which of the vacua is the true vacuum. 
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IRX 



IRU 



UVU 




12.4 



uvx 




UVU 



(a) 



(b) 



FIG. 1. Solution of the SDE at p/ f n = 0.3 for Njj = 70, Nx = 60: (a) real part; (b) imaginary 
part. Scale Aqcd 1S indicated by on U and X axes. Each number on the top of the vertical axis 
is the maximum value of HeB/ f n or ImB/f^. 



Integrand of Xf(M-) Integrand of -V 




(a) (b) 

FIG. 2. Integrands of (a) and (b) -F[-B so i] at p/fa = 0.3 for N v = 70, N x = 60. The 

upper 9/10 of each figure is clipped. 
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(40,30) o 

(50,40) + 

(60,50) □ 

(70,60) x 



FIG. 3. Typical values of fn{n)/fn for four choices of the size of descretization, 
(Nu,N x ) = (40,30), (50,40), (60,50) and (70,60). 



VZB S0 J#f 



(40,30) « 

50,40) + 

(60,50) □ 

(70,60) x 



FIG. 4. Typical values of V[B so {\/ f% for four choices of the size of descretization, 
(Nu,N x ) = (40,30), (50,40), (60,50) and (70,60). 
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FIG. 5. /x-dependence of /-n- (/i) /A at T = for two choices of the initial trial functions (A) 
and (B) in Eq. ( pi)) . 



As we discussed in Sect. O, the true vacuum is determined by studying the effective 
potential. We show the value of the effective potential V[B so i]/ f£ in Fig. || Since the value for 
the trivial solution is already subtracted from the expression in Eq. (|2.18|) , positive (negative) 
value of V^-Bgoi] implies that the energy of the chiral broken vacuum is larger (smaller) than 
that of the symmetric vacuum. For 3.0 < /i/A < 5.23, V^-Bsd]^) < V[B so \ 



Kb) 



0, which 



implies that the chiral broken vacuum is energetically favored. On the other hand, for 
5.23 < /VA < 6-8, V[5 S oi](A) > 0, and the chiral symmetric vacuum is the true vacuum. 
The chiral symmetry is restored at the point where the value of V[-B S oi] becomes positive: 
the chiral phase transition occurs at /i = 460 MeV (/i/ f n = 5.23). Since the value of the pion 
decay constant vanishes discontinuously at that point, the phase transition is clearly of first 
order. The non-trivial solutions for 5.23 < \ij < 6.8 correspond to meta-stable states, 
which were shown in Ref. Jl3[ by assuming a momentum dependence of the mass function. 



The result here obtained by solving the SDE (2.15) agrees qualitatively with their result. 



C. T + and fj, = 



At non-zero temperature we perform Matsubara frequency sum by truncating it at some 
finite number. Integration over the spatial momentum is done as shown in Eqs. ( |4.1| ) and 
( |4.2| ). We have checked that the integration range given in the first line of Eq. ( |4.8| ) is large 
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FIG. 6. /i-dependence of the effective potential V[B so i}/ f£ at T = 0. 

enough for the present purpose. Since the trial function (B) converges to the same solution 
as (A) for small fi in the case of T = 0, it may be enough to use the trial function (A) for 
studying phase structure. Here and henceforth we use the trial function of type (A) only. 
First, we check dependences of the results on the truncation. In Fig. [7] we show the values 
of f x (T)/f n at T/fn = 0.5, 1.0 and 1.5 for three choices of the truncation point: Ny = 10, 
15 and 20. It is clear that the choice Njj = 20 is large enough as the truncation point for 
the present purpose. 

Resultant temperature dependences of fn{T)/ f n and V{B so \]/ are shown in Fig. |8|. 
The value of f^iT) does not change in low temperature region. It once increases around 
T/f n ~ 0.7 and decreases to zero above that point, and finally reach to zero around T//V = 
1.89. Differently with the previous case the value of the potential reaches to zero around 
the temperature where the decay constant vanishes. Since there are numerical errors in 
the present analysis, we can not clearly show whether the decay constant and the potential 
vanish simultaneously. Our result shows that the chiral phase transition is of second order 
or of very weak first order, and that the critical temperature is around T c = 166 MeV. 



D.T^O and \i + 

Now, we solve the SDE when both T and /i are non-zero. Infrared and ultraviolet cutoffs 
for x-integration are fixed as in the first equation of Eq. ( |4.8|) , and the sizes of descretizations 
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(10,60) o 
(15,60) + 
(20,60) □ 



FIG. 7. Values of U(T)/f n at T/f w 
point: N v = 10, 15 and 20. 



0.5, 1.0 and 1.5 for three choices of the truncation 



are fixed to be (N U7 N x ) = (20, 60). 

In the previous subsections we found that the phase transition at T = and /i ^ is 
of first order, while that at T ^ and fi = is of second order or of very weak first order. 
Then one can expect phase transitions of first order for small T and large fi, and those 
of weak first order for large T and small /i. First, we show the temperature dependence 
of U{T)/U m Fig. | and that of V[B sol }/f^ in Fig. [TO| for = 1, 2, 3, 4 and 5. 

Figure |9| shows that in all the cases the values of the pion decay constant once increase 
around T ~ f n and decrease above that. These, especially for fi/f 7T = 1 and 2, behave as if 



the phase transitions are of second order. However, it is clear from Fig. |10| that the values of 
the effective potential V[-B so i] become positive before the values of the pion decay constant 
vanish. Then the phase transitions are clearly of first order as in the case of T = and 
7^ 0. To study phase transitions for smaller /i we concentrate on temperature dependences 
around phase transition points. Shown in Figs. |ll] and [L2| are temperature dependences of 
the pion decay constant and the effective potential for jjl/ f n = 0.25, 0.5 and 0.75 together 
with those for fj,/ f n = 1. For \ij f n = 0.5, 0.75 and 1.0 the values of the pion decay constant 
approach to zero as if the phase transitions are of second order. However, the values of the 
effective potential become positive before the values of the decay constant reaches to zero. 
Then we conclude that the phase transitions for /a/ f n = 0.5, 0.75 and 1.0 are clearly of first 
order. On the other hand, the phase transition for fi/f n = 0.25 is of second order or of very 
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f„m/f n 



© o o o 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



(a) 



(b) 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



FIG. 8. Temperature dependences of (a) f n (T)/f n and (b) V[B so i]/f%. 
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0.2 0.4 0.6 0.8 



1.2 1.4 1.( 



FIG. 9. Temperature dependences of J^l/ fir for fi/f^ = 1, 2, 3, 4 and 5. 




FIG. 10. Temperature dependence of the effective potential V[B so \]/f% for /j,/f w = 1, 2, 3, 4 
and 5. 
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1.76 1.78 
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0.75 + 

0.5 □ 

0.25 x 



FIG. 11. Temperature dependences of U(T)/f n for \xjU = 0.25, 0.5, 0.75 and 1.0. 



1.76 1.78 



1 o 

0.75 + 

0.5 □ 

0.25 x 



FIG. 12. Temperature dependence of the effective potential V[B so \]/ f% for fi/ f w = 0.25, 0.5, 
0.75 and 1.0. 
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weak first order. 



T (MeV) 

200 | 1 



50 100 150 200 250 300 350 400 450 500 

(MeV) 

FIG. 13. Phase diagram obtained by the present analysis. Points indicated by o are the phase 
transition points of first order, and points by + of second order or very weak first order. 



Finally we show a phase diagram derived by the present analysis in Fig. [U| As was 
expected, phase transitions for small T and large \i are of first order, and those for large T 
and small // are of second order or of very weak first order. 



V. SUMMARY AND DISCUSSION 



We analyzed the phase structure of QCD at finite temperature and density by solving 
the self-consistent Schwinger-Dyson equation for the quark propagator with the improved 
ladder approximation in the Landau gauge. A pion decay constant was calculated by using 
a generalized version of the Pagels-Stokar formula. Chiral phase transition point was deter- 
mined by an effective potential for the propagator. When we raised the chemical potential 
\i with T = 0, the value of the effective potential for the broken vacuum became bigger than 
that for the symmetric vacuum before the value of the pion decay constant vanished. The 
phase transition is clearly of first order. On the other hand, for T ^ and fi = the value 
of the effective potential for the broken vacuum reached to that for the symmetric vacuum 
around the temperature where the value of the pion decay constant vanished. The phase 
transition is of second order or very weak first order. We presented the resultant phase dia- 
gram on the general T — fj, plane. Phase transitions are clearly of first order in most cases, 
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and for small fi they are of second order or very weak first order. Our results show that it is 
important to use the effective potential to study the phase structure at finite temperature 
and density. 

Finally, some comments are in order. Generally the running coupling should include 
the term from vacuum polarization of quarks and gluons at finite temperature and density. 
Moreover, the gluons at high temperature acquire an electric screening mass of order gT. |28| 
We dropped these effects, and used a running coupling and a gluon propagator of the same 
forms at T = fi = 0. (The same approximation was used in Ref. fl7j .) These effects can be 
included by using different running coupling and gluon propagator which depend explicitly 
on T and /i such as in Refs. [jT~3| , pT6[ . 



The approximation of A — 1 = C = might not be good for high tempreture and/or 
density. However, inclusion of the deviations of A — 1 and C from zero requires the large 
number for truncating the Matsubara frequency, and it is not efficient to apply the method 
used in this paper. New method to solve the SDE may be needed. We expect that the 
inclusion of the deviations of A — 1 and C from zero does not change the structure of the 
phase transition shown in the present paper. 
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